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Gluonic profile of the static baryon at finite temperature. 
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The ghion flux distribution of a static three quark system has been revealed at finite temperature in 
the pure SU(3) Yang-Mills theory. An action density operator is correlated with three Polyakov loops 
representing the baryonic state at a temperatures near the end of the QCD plateau, T/Tc ~ 0.8, and 
another just before the deconfinement point, T/T^ « 0.9. The flux distributions at short distance 
separations between the quarks display an action-density profile consistent with a rounded filled 
A shape iso-surface. However the A shape action iso-surface distributions are found to persist 
even at large inter-quark separations. The action density distribution in the quark plane exhibits 
a nonuniform pattern for all quark separations considered. This result contrasts the well-known 
Y-shaped uniform action density gluonic-fiux profile obtained using the Wilson loop as a quark 
source operator at zero temperature. We systematically measure and compare the main aspects of 
the profile of the fiux distribution at the two considered temperature scales for three sets of isosceles 
triangle quark configurations. The radii, amplitudes and rate of change of the width of the fiux 
distribution are found to reverse their behavior as the temperature increases from the end of the 
QCD plateau towards the deconfinement point. Remarkably, we find the mean square width of the 
flux distribution shrinks and localizes for quark separations larger than 1.0 fm at T /Tc ~ 0.8 which 
results in an identiflable Y-shaped radius profile. Near the deconfinement point, the action-density 
delocalizes and the width broadens linearly at large quark separations. 

PACS numbers: 12.38. Gc 



I. INTRODUCTION 

Revealing the color field distribution in the nucleon is a 
subject of fundamental importance to quantum chromo- 
dynamics (QCD) and confinement. Lattice QCD simula- 
tions provide a first principle source of knowledge about 
how the energy distribution manifests itself among a sys- 
tem of three static quarks (3Q). This has to do with the 
relevant ansatz that accurately parametrizes and models 
the non-abelian force that binds the nucleon. The color 
distribution due to a 3Q system has been of a problem 
of reviving interest of lattice simulations and has been 
revisited with variety of lattice techniques [ll-Ql- Even 
though, an important aspect of this problem yet remains 
to be thoroughly investigated. That is, the energy dis- 
tribution associated with the 3Q system at finite tem- 
perature. Tackling the problem of the gluonic distribu- 
tion from this perspective involves the employment of a 
methodologically different set of unbiased hadronic oper- 
ators. In addition to that, revealing the changes of the 
gluonic profile of the (3Q) system under various tempera- 
ture conditions would certainly contribute to our percep- 
tion of the underlying gluonic picture and the associated 
gluon dynamics. In fact, the distribution of gluonic fields 
in the baryon at high temperature, before quantum chro- 
modynamics (QCD) undergoes a phase transition, is un- 
known in detail and has not yet been under the scrutiny 
of the lattice approach. 

Most of our current understanding of the (3Q) con- 



• ^abakry @physics . adelaide . edu. ail] 



fining force is based on the analysis at zero tempera- 
ture [5l4ll|. The parametrization which provides the 
best possible fits of the lattice data of the measured 3Q 
system potentials has been controversial for a long period 
of time [a, 0, [l2,[3- However, recent lattice QCD find- 
ings regarding the three quark potential are settled to 
support the so-called A-ansatz parametrization for small 
quark separation distances of _R < 0.7 fm and the Y- 
ansatz for 0.7 < R < 1.5 fm [l3|. The A ansatz accounts 
for a confining potential built up as a sum of two-body 
forces, the string tension is half that in the corresponding 
QQ system and the confining part of the baryonic poten- 
tial is in proportion to the perimeter of the triangle set 
up by the 3Q system. On the other hand, if the confining 
potential is proportional to the sum of the distances from 
the quarks to the Format point with a string tension the 
same as that in the QQ system, then due to its shape, 
this potential is known as the Y-ansatz, giving rise to a 
three-body term relevant to a genuine interaction channel 
of the non-abelian force. 

Ambiguities are known to arise, however, in the cal- 
culations of the gluonic distribution in the 3Q system at 
zero temperature. The energy distribution may be vul- 
nerable to systematic errors associated with excited-state 
contamination [2] when constructing the static baryon 
using a Wilson loop operator. The configuration of the 
spatial links that best minimize the potential has to be 
adopted before hand. Associated with this arbitrariness 
in tuning the ground state operator are the excited state 
potentials which manifest themselves in the revealed glu- 
onic profiles as a bias refiecting the form of spatial links 
of the Wilson loop operator ^. The L shape baryon 
operator provides a pronounced evidence where the flux 



distribution mimics the source Q. 

The isolation of the ground state is challenging in the 
case of field-distribution calculations which involve four- 
point correlations rather than the ground state potential 
which is extracted in the large time limit of a three-point 
correlation [2|- For example, Euclidean time evolution 
in the three quark Wilson loop operator results in ob- 
servable broadening of the junction in the Y-shaped con- 
figuration [3|. Statistical noise, nevertheless, imposes a 
practical constraint on any further increase in the expo- 
nentially decaying operator. 

In this investigation, the static baryonic states are ac- 
counted for by means of Polyakov loops. This provides a 
gauge invariant operator which acquires a methodologi- 
cal importance [J] due to the ability to construct an unbi- 
ased 3-quark operator without recoursing to a particular 
assumption regarding the form of the configuration of 
the spatial links in Wilson loops or the ultraviolet prop- 
erties of these parallel transporters [l^ . While carrying 
out energy density calculations into the zero temperature 
regime requires substantial numerical simulations with 
regard to the CPU time as well as the memory storage, 
the use of these stringless hadronic operators for reveal- 
ing the energy distribution at finite temperature is still 
an attractive idea from the practical feasibility point of 
view. This can be studied in conjunction with the ther- 
mal effects. 

At finite temperature, pure Yang-MiUs SU(3) lattice 
simulations for the action density in the mesonic sec- 
tor display a flux distribution with a vibrating string- 
like shape. The density distribution shows a non-uniform 
pattern with an almost constant cross section in the in- 
termediate distance region 0.5 < i? < 1 fm and non- 
constant cross section at larger quark separations [16|. 
The non-uniformity of the action density coincides with 
only a small decrease in the QQ effective string tension 
a [131 suggesting the ground state may also display a non 
uniform action density distribution [l5| . 

In this paper, we generalize this analysis to the dis- 
tribution of the color field inside the baryon. We con- 
sider one temperature near the end of the QCD plateau 
region at T /Tc ~ 0.8, and other just before the de- 
confinement point at T/Tc ~ 0.9. The three infinitely 
heavy quarks are accounted for by means of Polyakov 
loops of the same time orientation. The field strength 
inside the corresponding quark system is revealed by 
correlating an improved action density operator [18[ to 
these gauge-invariant hadronic operators. Gauge-field 
smoothing [l3|. in addition to a high statistics gauge- 
independent [3| averaging is employed to enhance the 
signal to the noise. This noise reduction approach can 
be employed in a controlled and systematic manner that 
has been proved effective in keeping the physics intact in 
the case of the static meson [1^] . The analysis on either 
the QQ force or the action density shows that smearing 
leaves no effect on the corresponding measurements taken 
for quark source separation distance scales greater than 
the diameter of the Brownian motion of a diffused link 



i.e. the characteristic diameter of smearing. Moreover, 
the systematic effects associated with this UV filtering 
procedure on the gluonic profile has been reported in de- 
tail in Ref. 0. 

The analysis on the smearing effects is revisited in this 
work for the 3Q force. The relevant distance scale where 
the physics is preserved is established. After identify- 
ing this scale, the characteristics of the action density 
profile are presented for selected 3Q configurations and 
contrasted at the two considered temperatures. 

The map of this paper is as follows: In Section HIl and 
Section lllll the details of the simulations and noise reduc- 
tion techniques are described. The force in the 3Q system 
for selected configurations is evaluated in Section ITVl In 
Section |Vl the main aspects of the gluonic profile of the 
baryonic action density is analyzed and contrasted at the 
two temperatures. In the last Section IVIl the conclusions 
are provided. 



II. MEASUREMENTS 

The infinitely heavy quark state is constructed by 
means of Polyakov loop correlators. In the confinement 
phase, for pure SU(3) gauge configurations, the correla- 
tors respect the center symmetry transformation 



Ui_i.=i{x,nr = 1) = CUi_,=i{x,nr = 1), 



(2.1) 



where center C of the gauge group 51/(3) is all the el- 
ements z such that zgz~^ = g, with g E SU{3) or 
z = cxp{2Tril/3) G Z{3) with I = 0,1,2. The form of 
the center symmetry preserving baryonic correlators is 
then 

(P3Q(ri,r2,f3)) ^ (P(r-i)P(f2)P(f3)) 

= (e2-'P(ri)P(r2)P(n,)> 
= (P(ri)P(r2)P(r3)) , 

where the Polyakov loop is given by 



P(r) = -Tr 
^ ' 3 
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(2.2) 



which corresponds to three Polyakov lines all in the same 
time direction. 

After the construction of the gauge-invariant color- 
averaged quark states, subsequent measurement by a 
gauge-invariant action density operator S{p,t) is taken 
at the spatial coordinate p of the three dimensional torus 
corresponding to each Euclidean time slice. The ac- 
tion density operator is calculated via a highly-improved 
©(af) three-loop improved lattice field strength ten- 
sor [13. 

A scalar field that characterizes the gluonic field can be 
defined as 



C{p,ri,r2,r3) 



{'P3Q{ri,r2,r3)S{p)) 

(7'3Q(ri,r2,f3))(5(/?)) 



(2.3) 



for baryonic systems, where < > denotes averaging 

over configurations and lattice symmetries, the vectors f^ 
define the positions of the quarks and p the position of 
the flux probe. Cluster decomposition implies C — > \ 
away from the quarks. 

In this investigation, we have taken 10, 000 measure- 
ments at temperature T /Tc = 0.8, and 6, 000 measure- 
ments at temperature T/Tc — 0.9. The measurements 
are taken on hierarchically generated configurations. The 
gauge configurations are generated using the standard 
Wilson gauge action on lattices with a spatial volume of 
36^. Gauge configurations are generated with a coupling 
value of /3 = 6.00. The lattice spacing at this coupling 
is a = O.I fm [2i|. After each 1000 of updating sweeps, 
?^sub ~ 20 or 12 measurements separated by 70 sweeps of 
updates are taken for the two lattices corresponding to 
T/Tc w 0.8 and T/Tc « 0.9 respectively. These sub mea- 
surements are binned together in evaluating Eq. (j2.3p . 
The total measurements are taken on 500 bins. 

The gluonic gauge configurations are generated with a 
pseudo-heatbath algorithm [22]. The heatbath is imple- 
mented by (FHKP) [2j, [2J] updating on the correspond- 
ing three SU(2) subgroups. Each update step consists of 
one heatbath sweep and 4 micro-canonical reflections. 



III. STATISTICS 

Gauge-independent noise reduction can be employed 
by making use of the space-time translational invariance 
of the hypcr-toroid. By space-translational invariance, 
measurements of the action density operator is repeated 
for each translated Polyakov loops correlator. The mea- 
surements are then averaged over the spatial volume of 
the lattice. The periodicity in the time direction also 
allows averaging over the time direction 
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(3.1) 



where Nt is the number of time slices of the lattice. The 
symmetry of the quark positions can be also exploited to 
gain a reduction in the statistical uncertainties. The flux- 
distribution has been averaged around all the symmetry 
planes of a given quark configuration. 

Local action reduction by smearing the gauge links has 
been performed on the whole four dimensional lattice. 
This procedure can be applied for correlating operators 
with polyakov loops. For example, the correlations with 
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the topological charge has been studied in Ref. 
ing the Cabbibo-Marinari cooling. Smearing the gauge 
field can be particularly helpful in reducing the statistical 
noise associated with evaluating the Polyakov loop cor- 
relators. However this step may result in the elimination 
of short distance physics and one has to be careful with 
regard to the number of smearing sweeps and the rel- 
evant distance scale where the physical observables are 
extracted. In the next section, we study the effects of 



gauge field smoothing on the 3Q force and determine the 
distance scale and corresponding smearing level where 
the physical observables are left intact. We smooth the 
gauge field with an over-improved stout-link smearing al- 
gorithm |26j . The scheme of the over improvement of this 
algorithm is such that a minimal effect on the topology 
of the gauge field [20| is ensured. 

In the standard stout-link smearing [13|, all the links 
are simultaneously updated. Each sweep consists of a 
replacement of all the links by the smeared links. 



i/p(a;) ^ exp(i(3^j(a;)) U^ix) , 



(3.2) 



with 
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(nlix) ^ n^ix)) 

tr(f]t(x)-fi^(x)). 



and 



n^ix)^ K]p^,I]i,,(a;) [/t(x), 
\'^^^^ J 

where I]^,y (x) denotes the sum of the two staples touch- 
ing U^{x) which reside in the p, — v plane. In the over- 
improved algorithm, however, the square staple is re- 
placed by a combination of plaquette and rectangular sta- 
ples. This ratio is tuned by the parameter e [26[. In the 
following we use a value of e = —0.25, with p^ = p = 0.06. 

Despite the noise reduction at short distance quark 
separation, the signal drawn in statistical noise at large 
distances if one opt merely to the link integration 
method [23 in evaluating the Polyakov loop s in Eq. (|2.3p . 

While the link integration method [23] provides effi- 
cient noise reduction at small quark separations for the 
expectation values of Polyakov loops in Eq. (|2.3p . noise 
remains problematic at large separations. The method 
is expected to be efficient for the QQ source separation 
range i? < 0.5 [2^ but has to be supplemented with a 
large number of measurements for the calculations of the 
larger distance QQ potential [l^l • The number of mea- 
surements has to be increased significantly for the cor- 
responding 3Q potential calculations with a three point 
correlator. The action density calculations would be, on 
the other hand, even more challenging in both cases j20| . 
The exponential reduction of noise provided by the lev- 
eling approach of the Liischer Weiss (LW) algorithm re- 
quires the hierarchical integration to be carried out over 
time sub-slices larger than the dcconfinement temporal 
extent tgub > l/^c H- For the temperatures consid- 
ered here, the division of measurements into binned sub- 
measurements resembles a one level implementation of 
the LW method, with updating the last time slice. 



IV. FORCES IN THE STATIC BARYON 

Unlike the force which is extracted from a three point 
correlator, the flux characterization Eq. (|2.3p . involves 



a four-point correlation function in the numerator pre- 
senting additional challenges with respect to the signal 
to noise level. The lattice space-time and configuration 
space symmetries can be auxiliary in enhancing the signal 
to noise ratio, however, a four dimensional gauge smooth- 
ing has to be employed to obtain a signal. 

The force in the 3Q system is a physical observable of 
direct relevance to the properties of the underlying en- 
ergy distribution. In the following we consider the effects 
of the gauge smoothing procedure on the force experi- 
enced by a test color charge. This can give indications 
on the relationship between source separation distance 
and the number of smearing sweeps where the changes in 
physics is minimal. Similar techniques have been adopted 
in Ref. [SO] in the determination of the large distance QQ 
force in vacuum with different levels of hyperbolic (HYP) 
smearing [la|. In the following, we consider the evalua- 
tion of the force via three Polyakov loop correlators. 
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FIG. 1: Schematic diagram for the isosceles configuration of 
the 3Q system. Tlie large balls represent the motion of the 
diffused field of characteristic smearing radius of Ra centered 
at the quarks (small spheres). 



For several levels of smearing corresponding to Usw = 
{20,40,60,80}, we numerically evaluate the force on a 
test color charge, assuming the transfer matrix interpre- 
tation is preserved as justified in Ref. [20|, the 3Q Poten- 
tial can be identified via a three loop correlator as 

(^3q) = {V{ri)V{r2)V{rz)) 

= CXp(-V3Q(fi,f2,f3)/T). 

The force on the third quark Q3 for the isosceles triangle 
configuration illustrated in Fig. [1] is measured through 
the definition of the derivative on the lattice [3ll ^ 



Pqs - - 



dV{R-A) 



2aT 



dR 

f {P{0,0)P{0,A)P{R,A/2)) \ 



log 



y{P{0)PiO,A)PiR+l,A/2))J 



(4.1) 
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FIG. 2: The force for the isosceles 3Q configurations with 
base lengths (a) A = 0.6 fm, (b) A = 0.8 fm and (c) A = 1.0 
fm, respectively. The x-axis denotes the position R of the 
third quark. Smearing effects are manifest for R < 0.95 fm, 
R < 0.85 fm and R < 0.75 fm for A = 0.6 fm, A = 0.8 fm, 
and A = 1.0 fm. Only subtle smearing effects remain beyond 
these distance scales. 



TABLE I: The characteristic radii Bf = a/-R|, + ^42/4 for the 
baryonic system at each smearing level for each configuration 
of Fig. [3 



fused field characteristic radius is defined as 



Config 


A = 0.6 
Lf = 0.35 


A = 0.8 
Lf = 0.46 


A = 1.0 
Lf = 0.58 


TT-sw 


Rf Bf 


Rf Bf 


Rf Bf 


40 
60 
80 


0.65 0.63 
0.75 0.80 
0.85 0.90 


0.55 0.68 
0.65 0.76 
0.75 0.85 


0.45 0.67 
0.55 0.74 
0.65 0.82 



The numerical values of the force measured on smeared 
configurations are reported in Fig. [2] for three isosce- 
les bases, A = 0.6, 0.8 and 1.0 fm. The repeated mea- 
surements on the data sets corresponding to increasing 
smearing levels indicate, in general, invariance of the 
force experienced by the test charge Qs under smear- 
ing at large distances. The loss of short distance physics 
is pronounced at small values of R which decreases as 
we increase the length of the isosceles base quark config- 
uration. In the following, our consideration of different 
isosceles 3Q configurations enable a systematic identifi- 
cation of the distance scale beyond which a given level of 
smearing has little effect on the physical observables. 

Define i?i?(nsw) to be the minimal distance beyond 
which a smearing sweeps up to risw does not affect the 
force Eq. (|4.ip . The values of Rp can be read from Fig. [5] 
Table |T] summarizes the values of Rp for each isosce- 
les configuration and smearing level. Since the effects 
of smearing relate also to the length of the isosceles base, 
we list for comparison the values of the corresponding 
effective range, Bp, defined as Bp = \/Rp + A'^/4, and 
also the distance from a quark at the base of the triangle 
to the Fermat point of the triangle Lp ( see Fig. [Ij. 

Clearly the range Rp is decreasing with the increase of 
the length of the isosceles base. A, indicating that smear- 
ing around the charges residing on the isosceles base de- 
creases the force exerted on the color test charge Qs as 
the charges Qi and Q2 become closer as in Fig. [TJ In- 
spection of the corresponding values of the above defined 
Bp, on the other hand, show that the decrease of Rp 
with the increase of A is such that the length of Bp 
is approximately constant. To gain an insight to what 
these observations may imply, we study the character- 
istics of the Brownian motion of the diffused field, and 
also the analogous values of smearing threshold Rp{nsw) 
measured for the QQ system. 

The diffuse field is Gaussian distributed Q through a 
ball centered at position r whose evolution with a smear- 
ing time T, in a four-dimensional smearing scheme |16| is 
given by 



G(r;T) 



1 



(47rDT)2 



exp 



4:DT 



(4.2) 
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Jd^rG{r;Ty 
/d3rG(r;T) 



1/2 



^a^pcusw 



(4.3) 

The proportionality constant c scales the number of 
smearing sweeps Ug.-^] in the improved stout- link smearing 
algorithm defined above with respect to APE smearing as 
defined for instance in Refs. [3|, 115[. The calibration pro- 
ceeds via comparing the respective number of smearing 
sweeps in each smearing scheme with respect to a given 
threshold [SJl (the reconstructed action-density [l^l nor- 
malized to a single instanton action S/Sq). This yields 
a value of c = 6.15(3) [1Q|. With p = 0.06, the number 
of smearing sweeps in the improved stout-link smearing 
algorithm scales as half the number of the correspond- 
ing smearing sweeps in APE smearing with the smearing 
parameter a = 0.7. 

After identifying this characteristic smearing range, 
the values of i?s (n-sw) are compared to the corresponding 
values of Rp{nsw) for the QQ system [l^ in Table HIl 
Inspection of the values reported for the mesonic force 
unveils that Rp is roughly equivalent to twice the smear- 
ing radius Rg. This suggests that the mesonic force is 
invariant under the smearing operation applied on the 
whole four dimensional lattice as long as the fuzzed balls 
centered at the quark source links are non-overlapping. 
Similar analysis on the action density shows that the re- 
gion free of smearing effects obeys the same invariance 
criterion [l6|. 

The distance Bp describes the minimal distance from 
the quarks Qi,2 to Q3 for which the measured force is 
invariant under a given number of smearing sweeps. The 
values of Bp in TableUcompare favorably with the values 
calculated for the quark-antiquark system in Table |lll 
This indicates that the smearing effects are immaterial 
as long as the length of the isosceles side is such that 
the fuzzed balls around any of the color charges Qi,2 and 
that around the test charge Q3 are non-overlapping. 



TABLE II: The characteristic radii Ra and Rf at each smear- 
ing level for mesonic systems beyond which the force is unaf- 
fected by smearing [la ]. 



No 


.Sweeps _RF(fm) 


Rsiim) 


2Rsiim) 




40 
60 
80 


0.65 
0.75 
0.95 


0.38 
0.47 
0.54 


0.76 
0.94 
1.04 



with D describing the diffuseness of the field. The dif- 



For the smearing radii considered here, a slight overlap 
of the fuzzed ball around each quark Qi,2 on the base of 
the isosceles is seen to have no observable effect on the 
force experienced by the test charge Q3. This observa- 
tion does not exclude the possibility of the three body 
channel of the interaction in the (3Q) system. The locus 
of the center of interaction may still be outside the two 
overlapping spheres in the base of the triangle. 



In Tabic HI the distance from the quarks Q1.2 to the 
Fermat point is indicated for each configuration. Simple 
variation calculus shows that for an isosceles triangle, 
the position of the Fermat point does not depend on the 
height of the triangle, i?, and the locus is fixed merely 
by the length of the base of the triangle such that, xf = 
A/ {2^). 

In summary, for the quark position geometry consid- 
ered in this work, the above analysis on the measured 
values of the force among a system of three quarks for 
each smeared data set of configurations suggest a conser- 
vative distance scale 7 beyond which the confining force 
on a given source is unchanged to be 




7 



2i?,, 



(4.4) 



(a) R = 0.6 fm, rts- 



60. 



This restricts the number of smearing sweeps to be such 
that the characteristic diameter of smearing does not ex- 
ceed the distance between at most two quarks. To take 
into account the distance between other sources, an addi- 
tional conservative measure will be to keep the distance 
to the Fermat point from any of the quarks of a given 
configuration outside the radius of smearing. However, 
in the present case this may be immaterial since the dif- 
ferences in the force measurements are well within the 
statistical errors. Equation (|4.4p indicates the distance 
scales where a specific characteristic of the action density 
might be affected by gauge field smoothing. The corre- 
sponding effects of gauge field smoothing on the revealed 
gluonic profile will be discussed also on several occasions 
below. 




(b) _R = 0.8 fm, ns 



:60. 



V. ACTION DENSITY 

A. Flux Iso-surface Profile 

The flux strength is measured as the correlation be- 
tween the vacuum action-density, S{p,t), and a gauge- 
invariant operator representing the quark states as pro- 
vided by Eq. (j2.3|) . The action density operator 
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/^E^Tr(F^^r)2, 



(5.1) 



is calculated through the 0{a'^) improved lattice version 
of the continuum field-strength tensor [l^l 




(c) R = l.O fm, ris 



pimp 



E 



W, 






(5.2) 



^(i,i 



where Cui^ are ixi link products in the v^x plane and Wi 
are coefficients selected to remove 0{a^) and 0{a^) er- 
rors. The evaluation of this operator on smeared config- 
urations filters out the UV divergences around the quark 
positions. The origin of the co-ordinate system is placed 
at the middle between the two quarks (52,3 on the y axis 



FIG. 3: Surface plot (Inverted) of the flux distribution C{p) of 
Eq. p.3|l evaluated in the plane of the (3Q) system p{x, y, 0), 
for isosceles configuration of base length A = 0.4 fm and sep- 
aration distances (a) R — 0.6 fm, (b) R — 0.8 fm and 
(c) R = 1.0 fm, at T = 0.8 Tc. The spheres refer to the 
positions of the quarks. 




FIG. 4: The flux action iso-surface at the quark positions, 
plotted together with a surface plot for the density distribu- 
tion C{p), in the 3Q plane at temperature T = 0.9 Tc, for 
equilateral triangular configuration R — 1.1 fm and A = 1.0 
fm. 




(a) _R = 0.4 fm, n^ 



:40. 



(b) R = 0.5 fm, ns 



40. 




(c) R = 0.6 fm, ns- 



60. 



(d) R = 0.7 fm, ns. 



60. 



at positions p{0, ±^, 0) and at distance R from the third 
quark, Qs, at p{x = i?, 0,0). The quarks reside on the 
plane p{x,y, 0). 

On calculating Eq. ()2.3j) . we find C{p) < 1, and C ~ 1 
away from the quark positions. The density distribu- 
tion in the plane of the quarks is plotted in Fig. [3] In 
general, the action density distribution is non-uniformly 
distributed as revealed in Fig. [3] through \5\ The dis- 
tribution C{p{x,y, z = 0)) has an action density maxi- 
mal curve along the middle line p{x,y = 0,z = 0) be- 
tween the two quarks Qi,2- With the increase of source 
Qa separation, the peak point along the maximal curve 
C{p{x,y = 0,z = 0)) shows only subtle movement, re- 
maining near the Fermat point of the triangle. These 
results contrast the Wilson loop results at large separa- 
tions |3J where the action density assumes a constant am- 
plitude along each arm of the Y-shaped profile. A convex 
curvature in the contour plot of flux density is manifest in 
Fig. \Mjo). This also contrasts the density plots obtained 
using the Wilson loop where the flux density assumes a 
concave curvature. 

Figure |3] discloses the flux surface plot of C in the 3Q 
plane and associated iso-surface for an isosceles configu- 
ration corresponding to a base A = 1.0 fm at the temper- 
ature T/Tc = 0.9. The flux iso-surface displays a clear 
filled A shape distribution. By moving the third quark 
further away, i.e. by increasing R, the A shape is found 
to persist. The sequence of frames in Fig. [S] displays sim- 
ilar results, this time at T = 0.8 Tc. It is important to 
note that this geometrical form of the density plot man- 
ifests itself at a temperature near the end of the plateau 
of the QCD-phase diagram [34| where the string tension 




(e) R = 0.8 fm, naw 



:60. 



(f) R = 0.9 fm, ns 



= 80. 




(g) R= 1.0 fm, risw =80 



(h) R = 1.1 fm, nsw = 80 



FIG. 5: Surface plot in the plane of the 3Q system p{x, y,z = 
0) and iso-surface of the flux distribution C{p;fi,f2,f3) for 
the isosceles configuration with ^4 = 1 fm and the third quark 
separation distance R as indicated. T = 0.8 Tc. 
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(a) T = 0.8Tc, 



(b) T = 0.9Tc. 



FIG. 6: Compares the flux contour lines of the density distribution in the 3Q plane C for triangular base A = 1.0 fm and third 
quark separation R = 0.9 fm at (a)r = 0.8 T^ and (b) T = 0.9 Tc, in 2 = plane. 



has been reported to decrease only by a value around 

10% [13. 

The contour lines and iso-surface of the flux do not 
exhibit a significant change with the temperature scale. 
Similarly, the effects of smearing do not cause defor- 
mations of the iso-surface profile outside of smoothing 
the interpolation of the fiux-lincs. This has been ob- 
served [49[ at relatively short distances employing the 
link integration for evaluating Polyakov lines in the fiux 
strength characterization (|2.3p . Even though the num- 
ber of smearing sweeps selected for each graph are larger 
than the limits set by the invariance of the force analysis 
of the last sections, the rendered graphs are not sensitive 
to increased levels of smearing. 

The flux distribution acquires a non-trivial transverse 
structure along the lines perpendicular to the x-axis. 
The fits of the transverse distribution along the lines 
p{xi,y,0) is returning good x^ for a Gaussian distribu- 
tion with varying amplitudes and widths from the third 
quark Qa position, R, to the y-axis as shown for instance 
in Figs. |3] or [51 This symmetry about the j/-axis in the 
X — y plane also exists in the perpendicular z direction. 

In a mesonic system, the width of the flux distribution 
is cylindrically symmetric around the line joining the two 
quarks. However, the existence of a third quark away 
from the y-axis breaks the symmetry of the width profile 
and the measured widths perpendicular to the plane of 
the quarks do indeed differ from the widths in the quark 
plane. In the forth coming section, we focus on dissect- 
ing the profile properties of the flux distribution within 
the quark plane, while the asymmetry aspect ratios are 
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FIG. 7: The density distribution C'(p) for the isosceles con- 
figuration with the base, A — 1.0 fm, and height R = 0.8 
fm at T/Tc = 0.8 (risw ~ 60 sweeps). Data are plotted for 
the transverse planes a: = 1 to 2; = 6. The lines correspond 
to the Gaussian fits to the density in each plane p{xi,y,0). 
The highest amplitude lies close to the Fermat point plane 
X — 2.88 of this 3Q configuration. 



reported separately in the last section. 



TABLE III: The amplitude, Hy{xi) (scaled by a factor of 10^), of the flux distribution at each consecutive transverse plane Xi 
from the quarks forming the base, A, of isoscoles triangle. The measurements for base source separation distance A — 0.6, 0.8 
and 1.0 fm for the temperature T/Tc = 0.8 are indicated as a function of the third quark position, Q3. 



Plane 
Qs = R/a 



6 



8 



9 I = 10 X = 12 I = 13 



A=0.6 fm 

07 

08 

09 

10 

11 

12 

13 
A=0.8 fm 

07 

08 

09 

10 

11 

12 

13 
=1.0 fm 

07 

08 

09 

10 

11 

12 

13 



A 



4.12(2) 4.34(1) 4.28(1) 3.97(2) 3.48(2) 2.84(2) 

4.24(2) 4.53(2) 4.56(2) 4.36(3) 4.00(3) 3.49(3) 2.85(3) 

4.31(2) 4.63(3) 4.73(3) 4.63(4) 4.37(4) 3.99(4) 3.49(4) 2.85(3) 

4.35(3) 4.67(5) 4.82(5) 4.79(6) 4.62(6) 4.33(5) 3.95(4) 3.46(3) 2.83(2) 

4.36(5) 4.66(7) 4.84(7) 4.85(8) 4.74(8) 4.53(7) 4.24(5) 3.86(4) 3.39(3) 2.80(3) 

4.35(6) 4.6(1) 4.7(1) 4.8(1) 4.8(1) 4.59(9) 4.35(7) 4.06(4) 3.72(4) 3.29(3) 2.73(3) 

4.30(7) 4.4(1) 4.6(1) 4.7(1) 4.7(1) 4.5(1) 4.3(1) 4.06(6) 3.81(4) 3.53(4) 3.14(4) 2.64(3) 

4.71(5) 4.92(3) 4.80(1) 4.40(2) 3.80(3) 3.07(3) 

4.79(5) 5.08(3) 5.06(3) 4.78(3) 4.31(4) 3.71(4) 3.01(3) 

4.83(5) 5.14(4) 5.22(4) 5.03(4) 4.68(5) 4.22(4) 3.66(4) 2.98(3) 

4.84(5) 5.13(5) 5.27(5) 5.18(5) 4.93(6) 4.58(5) 4.13(4) 3.60(4) 2.95(4) 

4.84(5) 5.06(5) 5.24(5) 5.23(6) 5.08(6) 4.80(5) 4.44(4) 4.02(5) 3.51(6) 2.88(6) 

4.84(5) 4.92(5) 5.12(5) 5.18(6) 5.12(6) 4.92(5) 4.61(3) 4.23(6) 3.8(1) 3.4(1) 2.75(9) 

4.84(5) 4.70(4) 4.88(3) 5.00(3) 5.04(5) 4.92(5) 4.64(2) 4.27(8) 3.9(1) 3.5(1) 3.1(1) 2.6(1) 

5.17(9) 5.43(6) 5.33(3) 4.86(1) 4.15(2) 3.31(2) 

5.22(9) 5.56(7) 5.58(5) 5.25(2) 4.69(1) 3.99(2) 3.19(3) 

5.2(1) 5.58(9) 5.72(7) 5.53(5) 5.11(3) 4.55(2) 3.89(4) 3.14(4) 

5.2(1) 5.5(1) 5.8(1) 5.71(9) 5.43(7) 4.99(6) 4.45(6) 3.82(7) 3.08(7) 

5.26(1) 5.3(1) 5.7(1) 5.8(1) 5.7(1) 5.3(1) 4.9(1) 4.3(1) 3.7(1) 3.0(1) 

5.4(1) 5.1(1) 5.5(1) 5.7(1) 5.8(1) 5.6(1) 5.1(1) 4.6(1) 4.0(1) 3.4(1) 

5.7(1) 4.7(1) 5.0(1) 5.4(1) 5.6(1) 5.7(1) 5.3(1) 4.6(1) 4.0(1) 3.4(1) 



2.7(1) 
2.9(1) 



2.4(1) 



TABLE IV: The squared width. 



Ty, in lattice units, of the flux distribution as in Table HIH 



plane 


x = 1 


x = 2 


x = 3 


x = A 


x = 5 


x = & 


x = 7 


x = 8 


x = 9 


x = 10 


x = 12 


X=13 


O3 = R/a 


























A=0.6 fm 


























07 


9.7(1) 


9.6(0) 


9.3(0) 


9.0(1) 


8.6(1) 


8.2(1) 














08 


10.0(1 


9.9(1) 


9.7(1) 


9.5(1) 


9.2(1) 


8.7(1) 


8.2(1) 












09 


10.2(1 


10.3(1) 


10.2(1) 


10.0(1) 


9.7(2) 


9.3(2) 


8.8(2) 


8.2(1) 










10 


10.6(1 


10.8(2) 


10.6(2) 


10.5(2) 


10.3(2) 


9.9(2) 


9.4(2) 


8.8(2) 


8.0(1) 








11 


10.9(2 


11.5(3) 


11.2(3) 


11.1(3) 


10.9(3) 


10.6(3) 


10.2(2) 


9.5(2) 


8.6(1) 


7.8(1) 






12 


11.3(2 


12.5(5) 


12.1(4) 


11.8(4) 


11.6(4) 


11.3(4) 


11.0(3) 


10.3(2) 


9.4(2) 


8.4(2) 


7.7(1) 




13 


11.7(3 


13.9(7) 


13.2(6) 


12.7(6) 


12.4(6) 


12.3(6) 


12.1(5) 


11.5(3) 


10.4(2) 


9.2(2) 


8.3(2) 


7.9(2) 


A=0.8 fm 


























07 


12.2(2 


11.9(1) 


11.3(0) 


10.8(1) 


10.2(1) 


9.5(2) 














08 


12.4(2 


12.2(1) 


11.7(1) 


11.2(1) 


10.6(1) 


10.0(2) 


9.2(2) 












09 


12.8(2 


12.6(1) 


12.1(1) 


11.6(1) 


11.1(2) 


10.5(2) 


9.7(2) 


8.8(2) 










10 


13.2(2 


13.1(2) 


12.5(2) 


11.9(2) 


11.4(2) 


10.9(2) 


10.2(2) 


9.3(2) 


8.2(2) 








11 


13.6(2 


13.6(2) 


12.8(2) 


12.1(2) 


11.5(2) 


11.0(2) 


10.4(1) 


9.6(2) 


8.5(3) 


7.5(3) 






12 


14.1(2 


14.0(2) 


12.9(2) 


12.1(2) 


11.3(2) 


10.8(2) 


10.3(1) 


9.6(2) 


8.6(4) 


7.6(4) 


6.8(4) 




13 


14.8(2 


14.2(2) 


12.7(1) 


11.4(1) 


10.6(2) 


10.1(1) 


9.8(1) 


9.3(3) 


8.5(5) 


7.6(6) 


6.9(6) 


6.4(6) 


A=1.0 fm 


























07 


15.4(4 


14.6(3) 


13.7(1) 


12.9(0) 


12.2(1) 


11.4(2) 














08 


15.7(4 


14.8(3) 


13.8(2) 


13.0(1) 


12.3(0) 


11.5(1) 


10.6(2) 












09 


16.1(5 


15.0(4) 


13.8(3) 


12.9(2) 


12.1(1) 


11.5(1) 


10.6(2) 


9.5(2) 










10 


16.6(5 


15.1(5) 


13.7(4) 


12.6(3) 


11.7(2) 


11.0(2) 


10.4(2) 


9.4(3) 


8.1(4) 








11 


17.0(6 


14.9(6) 


13.3(6) 


12.0(5) 


10.9(4) 


10.2(4) 


9.6(5) 


9.0(6) 


7.9(6) 


6.6(6) 






12 


17.2(7 


13.9(5) 


12.2(5) 


10.8(5) 


9.7(7) 


8.9(7) 


8.3(8) 


7.8(9) 


7(1) 


6(1) 


5(1) 




13 


17.1(9 


12(1) 


10(1) 


9(1) 


8(1) 


7(1) 


6(1) 


6(1) 


6(1) 


5(1) 


5(1) 


4(1) 
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TABLE V: The amplitude, Hy{xi) (scaled by a factor of 10^), of the flux distribution at each consecutive transverse plane Xi 
from the quarks forming the base, A, of isoscoles triangle. The measurements for base source separation distance A = 0.6, 0.8 
and 1.0 fm for the temperature T/Tc = 0.9 are indicated as a function of the third quark position, Q3. 



Plane 
Qs = R/a 
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9 I = 10 X = 12 I = 13 



A=0.6 fm 

07 

08 

09 

10 
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12 

13 
A=0.8 fm 

07 

08 

09 

10 

11 

12 

13 
=1.0 fm 

07 

08 

09 

10 

11 

12 

13 



A 



2.73(1 


2.83(1 


2.77(1 


2.58(1 


2.27(2 


1.88(2) 


2.71(1 


2.79(0 


2.75(1 


2.62(1 


2.42(2 


2.14(2) 


2.68(0 


2.73(0 


2.67(0 


2.56(0 


2.41(1 


2.24(1) 


2.64(0 


2.64(1 


2.55(1 


2.43(0 


2.30(0 


2.19(1) 


2.59(0 


2.55(2 


2.42(2 


2.27(1 


2.14(0 


2.05(0) 


2.53(1 


2.45(3 


2.29(3 


2.10(2 


1.95(2 


1.86(1) 


2.46(2 


2.35(3 


2.15(3 


1.93(3 


1.76(3 


1.66(3) 


2.88(2 


2.93(1 


2.85(0 


2.64(0 


2.33(1 


1.94(1) 


2.84(2 


2.86(2 


2.79(0 


2.63(0 


2.41(1 


2.14(1) 


2.80(3 


2.77(3 


2.68(1 


2.52(0 


2.35(0 


2.18(1) 


2.74(2 


2.66(4 


2.53(3 


2.36(1 


2.20(0 


2.07(0) 


2.67(2 


2.54(5 


2.38(4 


2.18(3 


2.00(2 


1.89(1) 


2.59(1 


2.44(5 


2.23(4 


2.00(4 


1.81(3 


1.69(2) 


2.51(1 


2.35(4 


2.11(4 


1.85(4 


1.64(4 


1.50(3) 


2.85(5 


2.84(4 


2.78(2 


2.59(0 


2.31(0 


1.96(0) 


2.79(6 


2.74(5 


2.68(3 


2.52(1 


2.33(0 


2.10(0) 


2.72(6 


2.61(6 


2.52(4 


2.36(2 


2.20(1 


2.06(0) 


2.64(6 


2.46(7 


2.33(5 


2.15(3 


1.99(2 


1.89(0) 


2.55(6 


2.31(7 


2.14(6 


1.93(5 


1.75(3 


1.66(1) 


2.46(5 


2.19(7 


1.98(6 


1.73(5 


1.54(4 


1.43(2) 


2.39(5 


2.12(7 


1.86(6 


1.59(5 


1.38(4 


1.25(3) 



1.79(2) 

2.01(2) 1.71(1) 

2.06(1) 1.89(1) 

1.99(1) 1.91(1) 

1.83(1) 1.81(0) 

1.63(2) 1.64(1) 



1.63(1) 

1.77(1) 1.56(1) 

1.77(1) 1.67(1) 1.49(1) 

1.64(0) 1.63(1) 1.56(2) 1.43(1) 



1.81(1) 

1.97(1) 1.70(1) 

1.96(1) 1.82(1) 1.60(1) 

1.83(0) 1.78(1) 1.69(1) 1.51(2) 

1.64(1) 1.64(0) 1.63(1) 1.58(2) 1.46(2) 

1.45(2) 1.46(1) 1.50(1) 1.52(2) 1.50(3) 1.42(3) 



1.82(1) 

1.91(0) 1.69(1) 

1.83(0) 1.74(1) 1.58(1) 

1.63(0) 1.63(0) 1.60(1) 1.48(1) 

1.41(1) 1.44(0) 1.49(1) 1.49(1) 1.43(2) 

1.21(2) 1.25(1) 1.32(1) 1.39(2) 1.43(3) 1.38(4) 



TABLE VI: The squared width. 



Ty, in lattice units, of the flux distribution as in Table HVl 



plane 
h = R/a 



X = 1 X = 2 X — S X — 4 X — 5 X — 6 x — 7 x — 8 x — 9 a; = 10a:: = 12a:: = 13 



A=0.6 fm 



07 


13.2(1) 


12.8(1 


12.5(0 


12.0(1) 


11.4(2 


10.9(3 














08 


13.5(1) 


13.0(1 


12.9(0 


12.6(1) 


12.1(1 


11.4(2 


10.7(3) 












09 


13.8(1) 


13.2(1 


13.2(0 


13.0(0) 


12.7(1 


12.1(2 


11.4(2) 


10.6(3) 










10 


14.2(1) 


13.3(2 


13.4(1 


13.4(0) 


13.2(0 


12.8(1 


12.2(2) 


11.3(2) 


10.5(3) 








11 


14.6(1) 


13.3(2 


13.5(2 


13.6(1) 


13.6(1 


13.4(1 


13.0(1) 


12.2(2) 


11.2(2) 


10.2(3) 






12 


15.1(2) 


13.3(3 


13.5(2 


13.7(2) 


13.9(2 


13.9(1 


13.7(1) 


13.1(1) 


12.1(2) 


10.9(2) 


9.8(2) 




13 


15.6(2) 


13.2(3 


13.4(3 


13.7(3) 


14.0(3 


14.3(3 


14.3(2) 


13.9(1) 


13.1(1) 


11.9(2) 


10.5(2) 


9.3(2) 


A=0.8 fm 


























07 


16.6(4) 


15.8(3 


15.1(1 


14.3(0) 


13.3(2 


12.4(3 














08 


16.9(4) 


16.1(3 


15.6(1 


14.9(0) 


14.0(1 


13.0(2 


12.0(3) 












09 


17.3(4) 


16.3(4 


15.9(2 


15.4(1) 


14.7(0 


13.8(1 


12.7(2) 


11.7(3) 










10 


17.9(5) 


16.5(5 


16.2(3 


15.8(2) 


15.4(1 


14.7(1 


13.7(2) 


12.5(2) 


11.3(3) 








11 


18.5(5) 


16.7(5 


16.4(4 


16.2(3) 


15.9(2 


15.5(1 


14.7(1) 


13.5(2) 


12.2(3) 


10.9(3) 






12 


19.3(5) 


16.8(6 


16.7(5 


16.6(4) 


16.5(3 


16.2(2 


15.7(2) 


14.7(1) 


13.4(2) 


11.8(3) 


10.5(3) 




13 


20.0(5) 


17.0(6 


17.0(5 


17.0(5) 


17.0(4 


16.9(4 


16.6(3) 


15.9(2) 


14.7(2) 


13.1(3) 


11.4(3) 


9.9(3) 


A=1.0 fm 


























07 


21.4(9) 


20.1(7 


18.8(3 


17.4(0) 


15.9(1 


14.5(3 














08 


21.9(9) 


20.4(8 


19.3(5 


18.1(2) 


16.7(0 


15.2(2 


13.8(3) 












09 


22(1) 


20.7(9 


19.7(6 


18.7(3) 


17.5(1 


16.1(1 


14.5(2) 


13.0(3) 










10 


23(1) 


21(1) 


20.2(7 


19.3(5) 


18.3(2 


17.1(1 


15.6(1) 


13.9(2) 


12.4(3) 








11 


24(1) 


22(1) 


20.7(9 


19.9(6) 


19.1(4 


18.1(2 


16.8(1) 


15.1(2) 


13.4(3) 


11.7(3) 






12 


25(1) 


22(1) 


21.4(1 


20.8(7) 


20.1(5 


19.2(4 


18.1(3) 


16.5(2) 


14.7(2) 


12.8(3) 


11.2(4) 




13 


27(1) 


23(1) 


22.4(1 


22.0(8) 


21.4(7 


20.7(5 


19.7(4) 


18.2(3) 


16.4(3) 


14.3(3) 


12.3(4) 


10.5(4) 
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B. Flux Radius profile 

At large quark separation, the revealed flux tube pro- 
file using the 3Q Wilson loop operator at zero temper- 
ature exhibits a uniform tube amplitude with a radius 
that is only slightly increasing up to the position of the 
junction K^. Although the bias of the revealed energy 
distribution by the shape of the configurations of the 
spatial links [2, 3 leaves these rendered energy distri- 
butions somewhat uncertain, this flux distribution has 
been considered consistent with the parametrization of 
the 3Q ground state potential with a F-ansatz at large 
distance Q- The Y-shaped gluonic distribution has also 
been considered in consonance with the dual supercon- 
ductivity picture J43l - l45| of the QCD vacuum. The flux 
is squeezed into a thin region dual to the Abrikosov vor- 
tex [4y| resulting in the formation of Y-shaped string-like 
flux tube di-Sl. 

At finite temperature, on the other hand, one intu- 
itively would expect the quantum vibrations of the under- 
lying three string system |36l438{ to give rise to a nonuni- 
form action density distribution in a similar fashion to the 
results revealed in the meson [1^. The thin string-like 
Y-shaped flux tube may delocalize away from its classical 
configuration and span the whole region though out the 
bulk of the triangular 3Q configuration, giving rise to a 
rounded concentric family of A action iso-surfaces (equi- 
action surfaces of Fig.|3|). Each surface is weighted by a 
temperature-dependent amplitude intensity distribution. 
In this non- uniform action density context, the radius 
topology is not fixed merely based on the distribution of 
Equi-action surfaces, as there can be an infinite number 
of iso-surface topologies of the action density that all cor- 
respond to the same measured square root of the second 
moment of the distribution. 

The second moment. ry{x), and the amplitude, Hy{x), 
of the fiux density at each line p{xi, y,0) is measured by 
means of Gaussian fits to the complementary distribution 
C = 1-C 



C'{p{x,,v,Q))=Hy{x,)e-y"/^-" 



(5.3) 



The fits to this Gaussian form are illustrated in Fig. [T] 
The mean square width in the 3Q plane at position, Xi, 
is measured via 



'rlixi) 



J dyy^C'ip{x,,y,0)) 
J dyC'{p{xi,y,0)) 



(5.4) 



eliminating dependence on , Hy{xi). The values of the 
measurements of . Hy{xi) and ry(xi) are listed in Ta- 
bles mil through IVIl The radius profile in the quark 
plane, 2 = 0, is measured at each lattice co-ordinate Xi. 
The data points corresponding to radii along the x-axis 
for a given quark configuration are interpolated with a 
continuous line up to the third quark, Q^, with position 
p(i?, 0, 0) as in Fig. E 

The first row of graphs in Fig. [5] correspond to radii 
measurements at the temperature T/Tc = 0.8 with base 



length running from A = 0.6 fm to A = 1.0 fm. For 
A = 0.6 fm, the radius profile draws almost constant lines 
with small declination indicating a subtle decrease along 
the X-axis up to the third quark position. The difference 
in radii between the very first planes and the planes close 
to the third quark Q3 becomes more pronounced with the 
increase of the third quark Qs separation i? as well as the 
increase of the distance between the two quarks Q\_i in 
the base. 

At the same temperature scale T jT^ = 0.8 and small 
isosceles base A = 0.6 fm, the tube's radius, ry{xi), at 
a given point broadens slowly with the increase of the 
quark separation Q^. This behavior changes as the length 
of the isosceles base A becomes wider. The change in ra- 
dius along the x-axis with R approaches near a stagnation 
in the broadening for A = 0.8 fm indicating an inflection 
point. 

This is evident from the profile at the widest base 
length A = 1.0 fm. The radius at a given point typically 
decays with the increase of the third quark Q3 separa- 
tion R. The geometrical area spanned by the triangle 
made up by the 3Q system becomes significantly large 
as the third quark Q3 steps farther away. In response, 
the gluonic energy condenses in narrower extents around 
the X-axis. In fact, we see a clearly identifiable Y-shaped 
profile of the gluon flux emerging at i? = 1.3 fm as well 
as at i? = 1.2 fm. 

For illustration, the action contours and radius pro- 
file have been superimposed in Fig. |9l The action con- 
tours are concentric convex A-shaped action isolines, the 
corresponding radii measured along the x-axis are not 
coinciding with any of the action isolines and have a con- 
vex Y-shapcd-like profile. Variation of the amplitude, 
Hy, hides the underlying Y-shape revealed through the 
consideration of Vy. We see from the corresponding last 
graph of the bottom row in Fig. |S] that the contour lines 
and the radius profile have similar concave curvatures at 
the temperature T/Tc = 0.9. 

This analysis may render the widely used terms such 
as the Y and A shaped gluon flux linguistically ambigu- 
ous if their usual usage is brought to the regime of non- 
uniform action density profiles with position-dependent 
(local) amplitudes distribution, Hy. One can speculate 
that the observation of the simultaneous co-existence of 
both the Y and the A aspects of the profile opens the pos- 
sibility that the ground-state baryon state may exhibit a 
similar action isosurface behavior even with the success of 
the Y-ansatz in the parametrization of the large distance 
potential. This may sound plausible especially if we take 
into account the fact that the gluonic junction broadens 
with the evolution of Euclidean time in the Wilson loop 
operator |3|. In addition, non- uniformly UV regulated 
Wilson loop operators have optimized the ground state 
overlap at zero-temperature in the mesonic sector jl5| . 
At the same time, and along the line of the above argu- 
ment, the observation of a Y-shape current distribution 
following abelian gauge fixing [J at T/Tc = 0.8 should 
not be taken by no means as contradictory to the A- 
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FIG. 8: The radius profile of tfie flux-tube measured in the plane of the quarks for each isosceles configuration with base 
A = 6a, A = 8a and A = 10a (a = 0.1 fm), at two temperature T/Tc = 0.8 (top) and T/To = 0.9 (bottom). The legend (in 
the upper right corner graph) signifies the third quark position. 




FIG. 9: The radius profile of the flux-tube displaying a Y- 
shape like profile for quark configurations of base A — 1.0 fm 
and the third quark position R — 1.3 fm at temperature 
T/Tc = 0.8. In the background is the corresponding flux 
action-density contours. 
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shaped action density in QCD observed without gauge 
fixing. 

Inspection of the bottom row of graphs in Fig. |8]reveals 
how thermal effects on the tube's radius profile take place 
as we get closer to the deconfincmcnt point at T/Tc = 0.9. 
In general, there is an increase in the radius of the flux 



FIG. 10: Same as Fig. [8] for isosceles configuration bases of 
A — 0.6 fm, and A — 1.0 fm. The top and bottom figures 
compare the measured radius profile for two levels of smear- 
ing, 60 sweeps and 80 sweeps, respectively. The radius in 
lattice units. 
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FIG. 11: The profile of the action density amplitude, Hy{xi) (scaled by a factor of 10^) for each isosceles configuration with 
base A = 0.6 fm, A = 0.8 fm and A = 1.0 fm, for the two temperatures T/Tc = 0.8 (top), T/Tc = 0.9 (bottom). The legend 
signifies the third quark position. 



with the increase of the temperature. The tube's topol- 
ogy is almost the same with an expansion of the size as 
we get to wider triangular bases. The radii flatten out 
through the planes a; < 6 for large quark separations. 
The change in radius along the x-axis increases also with 
the increase of the temperature for small quark separa- 
tions. Minimum growth in the radius for increasing R 
is noticeably manifesting near the Fermat point of the 
configurations x = 1.7, 2.3 and 2.9 for the isosceles bases 
A = 0.6,0.8 and 1.0 fm, respectively. Another distin- 
guishable feature for the profile at T/T,, = 0.9 is that the 
radius docs not show any sign of squeezing at all quark 
configurations. The increase in energy resulting from the 
increase of the temperature is now large enough to ac- 
commodate the corresponding enlargement in the geo- 
metrical area of the triangle set up by the quarks. We 
focus on detailed aspects of the flux broadening patterns 
separately in Sec. D. 

In addition to the force measurements in Sec. IIVI taken 
as a guiding analysis to set a trusted distance scale for 
each level of smearing, we now report the effects of smear- 
ing on the radius profile of the action density along the 
tube. Fig. [TUl compares the radii of the flux at each plane 
X measured on 60 and 80 sweep smeared gauge configura- 
tions. The values of the measured radii do not change at 
distant planes from the isosceles base. Smearing causes 
a subtle shifting rather than lensing effect on the radius 
at the planes near the quarks in the base Qi^2- An in- 
crement of 20 sweeps of smearing from 60 sweeps to 80 
sweeps causes a maximum increase of the radius by a sub- 
tle factor of 1.04. This effect diminishes as we consider 
far planes a; > 6 from the Qi^2 quarks on the base. 



C. Flux Amplitude profile 

At zero temperature, the revealed vacuum structure 
inside the static baryon constructed via Wilson loop op- 
erator has a maximum vacuum suppression at the center 
of the triangle made up by the 3Q system at small separa- 
tions 0]. At large distances, the Wilson loop operator of 
the minimum spatial string length has been found to min- 
imize the potential [3|, indicating a junction position at 
the Fermat point of the configuration. However, a peak 
in the action density at zero temperature does not man- 
ifest. The distribution assumes a constant amplitude. 
The analysis performed here for the density distributions 
using Polyakov lines as hadronic operators, nevertheless, 
reveals density amplitude peaks which manifest at short 
as well as large source separation distances. In this sec- 
tion, the amplitude profile is investigated and contrasted 
for the two considered temperatures. 

Tables Hill and FVl summarize the measured amplitudes, 
Hy{x), in accord with the Gaussian fits of Eq. (|5.3p . The 
corresponding plots are shown in Fig. Ill[ for the isosce- 
les configurations with base A = 0.6 fm, A = 0.8 fm 
and A = 1.0 fm at two temperatures T/Tc = 0.8 (top 
row) and T/Tc = 0.9 (bottom row), respectively. At all 
considered planes, the height of the distribution Hy(x) 
decreases with the increase of the temperature, which 
reciprocates the changes of the radius of the flux with 
the temperature. The decrease of the distribution height 
together with the associated increase in the distribution 
moment indicates the spread of the gluonic energy with 
the increase of the temperature. 

At T/Tc = 0.8, the amplitude also increases at most 
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planes when moving the third quark Q^ farther from the 
base of the isosceles configuration. Recalling the corre- 
sponding decrease in the radii along the x-axis, one infers 
the gluonic behavior undergoes a localization rather than 
a decay of the flux tube. The amplitudes at T /Tc = 0.9 
show similar increase up to small quark separation. How- 
ever, a noticeable turn over to decreasing amplitude with 
the increase of the third quark Qs separation manifests 
for R'> 10 a. The behavior of the amplitude and radius 
at T/Tc = 0.8 resembles respectively the behavior of the 
radius and amplitude at T/Tf, = 0.9. The analysis of the 
flux amplitude shows different qualitative behavior as we 
transit from the end of the QCD plateau to just before 
the deconfinement point and this behavior is reciprocal 
to the radius profile indicating a delocalization of the 
gluonic distribution with the increase of the temperature 
and the subsequent decrease in the string tension. 

In addition to the position of amplitude maximum 
along the x-axis at y = and the corresponding trigono- 
metric aspects of the triangular setup. The maxima lo- 
calize around the second and third planes x = 2,x — i, 
for third quark Q3, separations R < 10 a. However, this 
localization of the maxima of the vacuum suppression 
around the Fermat points ceases as the third quark is 
pulled away further. The density maximum moves in the 
same direction of third quark at T /T^ = 0.8 and moves 
in the opposite direction (towards the triangle base) for 
the higher temperature T/Tc = 0.9. 



D. The broadening of the flux width 

In this section, we focus on the broadening aspects of 
the mean-square width of the flux. We restrict our anal- 
ysis to the mean square width in the 3Q plane at the two 
considered temperature scales. The lattice data for the 
mean square width, r'^ixi), at planes Xi along the x-axis 
are summarized in Tables IIIII and |Vl For convenience, 
we have considered the tube's width for an analysis per- 
formed on gauge configurations of 80 sweeps of smearing, 
where we obtain the best signal to noise ratio with only a 
relatively small elimination of short distance points which 
are affected by smearing. 

In the last section, we reported the effects of gauge- 
smoothing on the radius of the gluon flux. Smearing 
shifts the width by a subtle constant near the base of 
the isosceles triangle, as in Fig. [121 This shift diminishes 
at distant planes from the base. To further examine the 
rate of broadening of the flux distribution, we fit the 
mean square width to the simple linear ansatz 



r;iR;x^) 



bi{xi)R + b2{xi) 



(5.5) 



The returned values of the slope of the growth in the fiux 
width with the increase of the isosceles height, R, display 
small systematic errors associated with the selection of 
parameters as indicated in Table I VIII 



TABLE VII: The slope of the growth in the mean square 
width, r^, measured for isosceles base A = 1.0 fm on two 
levels of link smearing. The measurements are obtained from 
the fits to the linear form Eq. (|5.5|l . 
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60 


0.76(4) 


4-13 


80 


0.83(3) 


4-13 


x=7 






60 


2.3(1) 


8-13 


60 


2.8(2) 


10-13 
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2.2(1) 


8-13 
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The profile of the broadening of the glue at various 
planes, for isosceles base, A = 0.8 fm, is plotted in Fig. [13] 
with a similar plot for, A ~ 1.0 fm, in Fig. [T31 respec- 
tively. Each set of data describes how the width of the 
gluonic fiux vary at a given plane Xi for the triangle base 
as the third quark Qa moves to a larger values of R. 
Evidently, the increase of the temperature dramatically 
increases the rate of the broadening of the glue at all 
planes. 

Apart from the pronounced thermal effects near the 
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FIG. 12: The squared flux distribution width at the depicted 
planes, a::^ = 2, 3 and 7, are compared for two smearing levels. 
The isosceles configuration base length A = 1.0 fm at tem- 
perature T /Tc = 0.9. Smearing merely shifts the profile by a 
constant. The broadening pattern is not afi^ected. 
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FIG. 13: The squared flux-tube width at the depicted planes for the isosceles configuration A = 0.8 fm compared at two 
temperatures T/Tc = 0.8 (top) and T/Tc = 0.9 (bottom). The plane coordinates are indicated in the legend. 
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FIG. 14: Same as Fig. [T3]for a larger isosceles base length of j4 = 1.0 fm. 



deconfinement point, we see the rates of broadening at 
T/Tc = 0.8 are decreasing as one proceeds to the more 
distant planes from the base of the triangle. Moreover, 
the wider the base of the isosceles triangle, the more pro- 
nounced is the corresponding decrease in the width, in- 
dicating that the gluonic tend to become more localized 



as the geometrical area enclosed by the quarks positions 
becomes larger. 

The shrinking of the width of the flux tube is a pe- 
culiar property to certain geometrical configurations of 
the Multi-quark system. The decrease in the width with 
the increase of the inter-quark separation has never been 



16 




© 



® 



T=0.8Tc, x=3 
A=6a ^•^ 
A=8a - o 



30 

25 

N 20 

CO 



10 11 12 13 



- 






T= 
A 


=6a ^ 


X: 


=3 








A: 


:8a ^ 


©- 


-' 


n 


® 





© 


$ 




^ o 


1 


• 


• 


• 


• 




• 




1 


1 


1 


1 




1 



10 
5 

30 
25 



■N 20 

CO 



10 
Ra-^ 



11 12 



13 



10 
5 



8 

8 



8 

9 



T=0.8Tc, x= 

A=6a ^-m- 
A=8a © 



S 
10 



1) 



11 



S 

12 13 



Ra-' 







T=0.9Tj,, x=7 
A=6a ^•^ 






A=8a - © 


i 


s 


o t °. - 
• • 









10 11 12 13 
Ra-^ 



FIG. 15: Comparison of the mean square width of the flux distribution at three distinct planes, a: = 2, 3, and 7, for two isosceles 
bases A = 0.6 fm and A = 0.8 fm. The upper graphs show the comparison at T/Tc — 0.9 whereas the lower are at T/Tc — 0.8. 



observed in the meson either using Polyakov lines at fi- 
nite temperature [l6| or Wilson's loop at zero tempera- 
ture [43] . The analysis of the Wilson loop based energy 
distribution at zero temperatures does not seem to indi- 
cate shrinking of the width of the flux tube 0] . 

Near the deconfinement point, the broadening of the 
mean square width, r^ , exhibits a clear linear divergence 
at distant planes from the base of the isosceles triangle 
for large separations, R > 1.0 fm of the third quark Q3. 
This result resembles the observed linear growing in the 
flux distribution width at the same temperature in the 
meson [20]. The slope of the growing width at distant 
planes from the base, a; > 6, from plane to plane show 
only subtle changes. This indicates that the effects of the 
boundary and the junction fluctuations decays away by 
proceeding to large quark separations. Fig. [15] displays 
the effects of increasing the distance between the two 
quarks (5i,2 on the rate of change of the width versus the 
motion of the source Q3. At T/Tc = 0.9 the wider the 
base the faster the rate of growth. This behavior is the 
reciprocal of the corresponding one at T/Tc = 0.8. 



E. Planes aspect ratio 

The gluonic flux in the 3Q system does not exhibit 
a symmetry between the width measured in the quark 
plane and that in the perpendicular direction. This has 
to do with the underlying gluonic structure and the as- 
sociated fluctuations. For example, the string picture 



indicates an asymmetry in the mean square width be- 
tween the two planes [301 • We report for completeness 
the general qualitative features of this ratio of the ac- 
tion density in the two perpendicular planes. The width 
of the tube in the perpendicular direction is measured 
through Gaussian flts as 



U^r) = 



J dzC'{p{xi,Q,z)) 



(5.6) 



We measure the ratio between the width in the quark 
plane and that in the perpendicular plane to the quarks 



a{xi) 



_ rl{x{) 



(5.7) 



It is interesting to consider this quantity since the pre- 
dictions of the string model at zero temperature indicate 
asymmetric width pattern near the junction [36|. In a 
similar way to the characteristics of the flux that we have 
studied in the previous sections, we plot in Fig. [16] the 
aspect ratio at the two temperatures for the same quark 
position conflgurations. 

Generally, the value of the aspect ratio indicates that 
the fluctuations in the quark plane are always larger 
than the perpendicular fluctuations for both tempera- 
tures. Further inspection of Fig. [TC] shows only subtle 
dependence on the temperature for the small isosceles 
bases. The aspect ratio is changing slowly as we move 
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to the quark plane for three isosceles bases A = 0.6 fm, A = 0.8 fm and A — 1.0 fm. The upper graphs show comparison at 
T/Tc = 0.8 whereas the lower at are T/Tc = 0.9. 



through the planes up to the thud quark position for 
the smallest isosceles bases. However at larger bases, the 
asymmetry throughout the gluonic cone becomes pro- 
nounced. The results of the aspect ratio are indicating 
greater restoring forces for the gluonic distribution in the 
perpendicular direction to the quark plane. This effect 
diminishes as we consider planes away from the Fermat 
point of the triangle quark configuration. Near the de- 
confinement point T/Tc = 0.9, the aspect ratio increases 
with the base length showing only a subtle dependency 
on R. This behavior is reversed at the lower tempera- 
ture T/Tc = 0.8 for large quark separation. Indeed we 
see at T/Tc = 0.8, for the bases A = 0.8 fm and A = 1.0 
fm, a decrease in the aspect ratio at all planes with the 
increase of the third quark separation i? > 1.0 fm. 

Within the underlying gluonic picture of three 
squeezed string-like flux tubes meeting at a junction and 
assuming a Y-shaped form, the predictions of this string 
model for the ratio of the mean square width of the 
flux distribution in the quark planes and the perpen- 
dicular direction at position of the junction has been 
worked out in Ref [36|. However, one should carefully 
consider the geometrical aspects of the configuration and 
take into account any remnant {T/Tc = 0.8) or pro- 
nounced {T/Tc = 0.9) thermal effects before pursuing a 
complete confrontation with the predictions of the bary- 
onic bosonic string models. We included the thermal 
effects [50| into the predictions that have been worked 
out earlier in Ref. [36[. However, this comparison of the 
lattice data with the predictions of the models lies be- 



yond the scope of the present presentation of the gluonic 
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profile, and will be reported in details elsewhere 

Finally, it is worth noting that the flux strength distri- 
bution revealed with the action density using the Wilson 
loop does not appear to produce an asymmetric gluonic 
pattern. For instance, in Ref. 3 the radius of the tube is 
calculated with cylindrical coordinates assuming a cylin- 
drical symmetry of the tube. The analysis provided here 
for the aspect ratio of the mean-square width provides an- 
other distinct feature of the glue, as revealed by Polyakov 
loops rather than a manifestation of temperature effects. 



VI. CONCLUSION 

The gluon flux distribution of a three quark system 
in pure SU(3) Yang-Mills vacuum has been revealed 
at finite temperature. This analysis is an extension 
of the calculations of the action density correlations 
obtained recently for the QQ [iQl system to three quark 
systems. The infinitely heavy baryonic state has been 
constructed by three Polyakov loops. The gluon flux is 
measured as a correlation between the action density 
operator and three traced (gauge-invariant) Polyakov 
lines. Measurements have been taken near the end of 
QCD phase diagram, T/Tc ~ 0.8 and just before the 
deconfinement point T/Tc ~ 0.9. 

Noise reduction has been achieved using a gauge- 
independent statistical approach that exploits the space- 
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time symmetries as well as symmetries of the quark con- 
figuration space. The calculations are performed at each 
point of the lattice and averaged over the lattice four- 
volume. The average over the configuration space have 
been calculated on 500 independent bins. Each bin con- 
sists of an average over measurements that are closely 
spaced in configuration space. An over-improved version 
of the stout-link smoothing algorithm has been employed 
with number of sweeps such that the physics is preserved 
in a systematic and controlled manner. 

We have revealed the characteristics of the fiux action- 
density measured for three sets of geometrical 3-quark 
configurations and the corresponding changes on the be- 
havior due to the temperature. Each set corresponds to 
isosceles triangle bases of length A = 0.6 fm, A = 0.8 fm 
and A — 1.0 fm. The characteristics of the isosurface, 
the radius and the amplitude profiles of the action den- 
sity correlations, in addition to the broadening (or the 
shrinking) pattern of the fiux distribution, can be sum- 
marized in the following main points: 

A. The Iso-surface of the fiux action-density dis- 
plays a family of concave A shapes at small as well as 
large quark separations. These A-shaped gluonic distri- 
butions persist and do not change into a Y-shape as the 
distances between the quark sources are increased. The 
density plots in the quark plane display a nonuniform 
distribution at all distance separations. This contrasts 
the Wilson loop results at zero temperature which 
exhibit uniform action density along each arm of the 
Y-shaped profile. A remarkable feature of the revealed 
map of the contour lines of the fiux strength is that 
the shape of the contour lines do not show significant 
sensitivity to the temperature for the two temperatures 
considered here. 

B. The Radius profiles give indications on the 
spread of the energy inside the baryon. At the lowest 
temperature near the end of the plateau, T/Tc = 0.8, 
the measurements of the radius indicate localization of 
the action density in narrow regions for quarks separa- 
tions greater than 1.0 fm. The radius of the tube de- 
creases and draws a Y-shaped like profile even though 
the action isosurface and isolines are A shaped. Near the 
deconfinement point, on the other hand, the energy tends 
to spread as we see the radius increases at all considered 
distance scales. 

C. The Amplitude profile analysis of the fiux den- 
sity shows a maximum vacuum fluctuations suppression 
at the plane nearest to the Format point of the pla- 
nar three-quark configurations for intermediate separa- 
tion distances. The distribution's peak ceases localiz- 
ing around the Fermat point of the 3Q isosceles con- 
figurations when the height, R, is greater than 1.0 fm. 
The peak shifts to the outside of the triangle made at 
T/Tc ~ 0.9 and shifts in the reverse direction to the in- 
side of the triangle for T/Tc ~ 0.8. That is the amplitude 
gets higher when the radius shrinks at T/Tc ~ 0.8 and 
the reverse manifest at T/Tc ~ 0.9. 



D. The Flux mean-square width does not always 
broaden with the increase of the quark source separation 
as is the case in the meson. For the lowest temperature, 
T/Tc ~ 0.8, the flux distribution shrinks in width for 
large quark separations. The change in the width of the 
flux tube shows a non-broadening aspect which is a prop- 
erty of certain configurations of the multi-quark system. 
The width, however, grows linearly near the deconfine- 
ment point, T/Tc ~ 0.9, with the increase of the height 
of the triangle. In general, the slope of the decrease or 
increase in the width, at both temperatures, depends on 
the length of the triangle base. The wider the base of 
the triangle set up by the quarks positions, the lower 
or higher is the slop at temperatures T/Tc = 0.8 and 
T/Tc = 0.9, respectively. 

E. The Aspect ratio between the mean square width 
of the fiux distribution in the quark plane and the width 
in the perpendicular plan exhibits an asymmetry. The 
gluonic fluctuations in the plane of the quarks are greater 
than that in the perpendicular directions around Fermat 
point, indicating a greater restoring force for the sys- 
tem in the plane of the quarks. The ratio between the 
two components of the mean square width decreases as 
we consider planes farther than the locus of the Fermat 
point of the quark configuration. The temperature de- 
pendence for the aspect ratio is more pronounced at large 
quark separations while we see almost the same profiles 
for small isosceles bases. The deviation of the aspect ra- 
tio from unity is implied by the predictions of the string 
models and do not manifest using Wilson loop operator 
in the action correlations. 

Future work 
As this work presents our first investigation of the flux 
distribution of the 3Q system at finite temperature, there 
are many promising avenues of investigation remaining. 
For example, additional quark configurations could be 
examined in detail. Lower temperatures remain to be 
investigated. We compare [S^ the growth of the mean 
square width of the flux distribution measured on the 
lattice with finite-temperature extensions of the baryonic 
string model [35j, |50| for the width of the baryonic junc- 
tion [3a]. On can also study the relevant ansatz for the 
measured potentials at each temperature. Methodologi- 
cal improvements based on increasing the number of mea- 
surements and decreasing the number of gauge smooth- 
ing sweeps are always desirable. The method pursued 
here may prove effective in a calculation framework that 
includes the effects of the dynamical quarks. 
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